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Abstract 



Under the Neumann constraints, each equation of the KdV hierarchy is decomposed 
into two finite dimensional systems, including the well-known Neumann model. Like in 
the case of the Bargmann constraint, the explicit Lax representations are deduced from 
the adjoint representation of the auxiliary spectral problem. It is shown that the Lax 
operator satisfies the r-matrix relation in the Dirac bracket. Thus, the integrabilitics 
of these resulting systems with the Neumann constraints are obtained. 

1 Introduction 

o 

CN ' It is well-known that quite a few finite dimensional integrable systems (FDISs) can be 

obtained from soliton equation by using the nonlinearization method of Lax pairs Q |j. 
Following [||], we call these FDISs the restricted flows, which can be divided into two kinds, 
<3\ ' one is the Bargmann system, which is a Hamiltonian system; and another is the Neumann 

system with some constraints ||. Many approaches have been developed to study the 
restricted flows. Among them is the r-matrix method. r-Matrix formula contains almost 
necessary information of the FDISs. In this paper, we develop the r-matrix method 
applicable for Neumann systems. We take restricted KdV flows as an example. We first 
^ review the results of the restricted KdV flows [Q, ||, || . A hierarchy of FDISs with Neumann 

constraints are obtained. We focus our attention on the FDISs related with KdV equation, 
including the well-known Neumann model H, which describes the motion of a particle on 
an (N — l)-sphere submitted to harmonic forces. We obtain their Lax representations 
from the adjoint representation of the auxiliary linear problem of KdV hierarchy. All the 
Lax matrices are 2x2 matrices (not N X N !). Then in Section 3, we show that the Lax 
operator L(X) satisfies the r-matrix relation in the Dirac bracket. Finally, by virtue of the 
general theory of r-matrix ||, we obtain the Uhlenbeck motion integrals Q and prove the 
involution property in the constrained space (sphere bundle). 
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2 The restricted KdV flows and their Lax representations 



It is well-known that the KdV hierarchy associates the following Schrodinger spectral 
problem 



■01 

i'2 



U(u,X) 



02 



U(X,u) 



1 

-X-u 



(1) 



To deduce the KdV hierarchy, following (8|, we first consider the adjoint representation 
of (S) 



V x = [U, V] = uv — vu 



with 



V 



( a ^ \ = y x~ m = sp ( a m b m \ 
V c —a I ' m ' \ c m —a m I 



which leads to 



b i+ i = Cbi, £ = -^d 2 - u + 1 u x , 



1 



bix j 



2oj+i — 2a>iU. 



The first few terms are 

a = b = 0, c = -1, a = 0, h = 1, ci 



«2 = ^i, ^2 = ~u, c 2 = ^ + u 2 ) , 6 3 = (u^ + 3u 2 ) , 
d 

where d = — — , a<9 _1 = c? _1 a" = 1. 
ax 

Then, we can take the following auxiliary spectral problem 



<M =V^(X,u) 



0i 

'02 



where 



n+l 



F(")(A,n) = VmX n+1 - m + 



m=0 







9n+2 







(2) 
(3) 

(4) 
(5) 

(6) 
(7) 



The compatible condition of (||) and (§) yields the zero-curvature equation 



u tn -v^+[uM n) 

which gives the KdV hierarchy 

Ut„ = — 26 n +2,x- 

For n = 0, = J7(«, A). 



0, 71 = 1,2,..., 



(8) 

(9) 

(10) 
(11) 
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For n = 1 we get 
/ 



V 



A 2 - -A + -4^ + — 



4 2 

and the KdV equation (set ii = i) 



1 

4 



(12) 



^(u xxx + 6uu 2 



(13) 



Now we turn to the nonlinearization of Lax pair for the KdV hierarchy. Following , 
we take N distinct X'jS and consider the following system 

ipij,x = ip2j, ip2j,x = (-Aj - u)ipij, j = 1, . . . , N. (14) 



It is easy to show that (up to a constant factor) 



5\i 



Ci/4j = AjV'ij. 



5u ^' 

As usual 0|, we can consider the constraints 
N x\ 

b k = y^, = 1,2, .... 



3=1 



In this paper we are interested in the Neumann constraint 
5u ' 



N X\ 

bl = j 



3=1 



i.e., 



1. 



Here and hereafter, we use the notation: qj = ipij, pj = ip2j, Q = (qi, 
(pi, . . . ,Pn) T and (•, •) denotes the stand inner product in space M. N . 
From (||) and (|l|) we get 

u = (p,p) - (Aq,q). 

Under this constraint, (|i~4|) is nonlinearized into the following system 



1x = P, 

p x = -Aq - ({p,p) - (Aq, q))q, 

iq) = h {q,p) = o. 

When N = 3, we have the case studied by Neumann 
call (|20| ) the Neumann model. 
Making use of the relation 



(15) 



(16) 



(17) 



(18) 



,Qn) T , P 



(19) 



(20) 

For convenience, we still 
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we have 

b i = (A i - 1 q,q), a i = -(A i - 1 q,p), c { = (A^W), i = 2,3,.... (21) 
Here and hereafter, symbol tilde means the corresponding nonlinearized quantity and 

/ N N 2 \ 

/ VhQh at \ 



V 



N 

£ 

k=l 



PkQk 
A — Xk 



V 



N o 

-i-E- ? * 



A' 

£ 

fc=i 
JV 

£ 



A — Xk 



PkQk 



(22) 



Under the condition fliDD, the following equation still holds 
V X = [U,V]. 



(23) 

Conversely a direct calculation can verify that ( p3| ) is just the Lax representation of equa- 
tion (|20|), For simplicity of notation, we set L(A) = —V, and Mj = Therefore we 
have 



Theorem 1. The equation (2C) has the Lax representation 
L x = [M ,L] 

with 

1 

-X+(Aq,q) -(p,p) 



L(X) 



A(X) 5(A) 
C(X) —A(X) 



M (A) 



where 



N 



m = £ 



Pkqu 
X- Xk 



N 



92 



A — Afc ' 



iV 2 

C(A) = 1 + E Pfc 



A — Afc 



(24) 
(25) 

(26) 



fe=i fc=i fc=i 

In the same way, we can discuss the nonlinearizations of the auxiliary spectral problem 



u 



j = l,...,N. 



(27) 



under the constraints (19) and (|2C|) and obtain a hierarchy of finite dimensional systems 
with the Neumann constraint (q, q) = 1. We call this the restricted t n -flow. Just like in 
the case of the Bargmann constraint, the identity 



V t , 



v {n \v 



(28) 



is a nonlinearization of the Lax representation of the restricted t n -flow. In particular, for 
n = 1 we obtain the following system 



1 1 

q t = Ap- (Aq,p)q+ - {p,p)p + -(Aq, q)p, 

p t = -A 2 q - ^(p,p)Aq + ^(Aq,q)Aq - (Ap,p)q 

+(^ 2 q,q)q - ^(A<?,<?) 2 <? + ^{p,p) 2 q + (&q,p)p, 

(q,q) = 1 



(29) 
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and the Lax representation 



Theorem 2. The equation (2i ) possesses the following Lax representation 



L X = [M U L], (30) 
M l( A)=| - {A<1 > P) X -\UP>P)- I (31) 



D(X) (Aq,p) 



where 



2 (p,p) - (Aq,q) , _\ , , A 2_ _\ (M,q) -(p,pY 



DW = -X 2 - " 1 A - (Ap,p) + <A^,g) 



Remark. After a simple calculation, we can show that: if (q,p) T satisfies the system of 



^^} = }^( 7 r- 7 r--^- 7 r- ■ (32) 



equation pQ ) and (29), then u = (p,p) — (Aq,q) solve the KdV equation (13). Besides 
from Q, we know it must be a finite-band solution of the KdV equation. Thus we can 
obtain the finite-band solution through solving the two finite dimensional systems with 
the Neumann constraints. 

3 r-Matrix for the restricted flow with Neumann constraints 

We know that the standard Poisson bracket for the two smooth functions f(p,q), g(p,q) 

N 

in the symplectic space M? N , ' S ^ j dpj A dqj is defined by 

N 

( 

j ^dqidpi dpidqi 
Since we are dealing with a constrained system with the holonomos constraints 

F= (q,q) -1 = 0, G = (q,p) = 0, (33) 
i.e., the constraint manifold is the sphere bundle TS N_1 

TS"- 1 = {(q,p) € R 2N \(q,q) = 1, (q,p) = o} (34) 

, we shall use the Dirac bracket ||], instead of the standard Poisson bracket. 
The Dirac bracket for the constraint (14) is || 

{/, 9}d = {/, 9} + \{f, F}{G, g} - ±{f, G}{F, g}. (35) 
A direct calculation yields 

{qk,qi}D = 0, {qk,Pi}D = hi ~ QkQh {Pk,Pi}D = -qkPi+Pkqi- (36) 



Proposition 1. The Neumann systems $2ui ) and ( pffi ) can be written as the following 
nonstandard Hamiltonian forms 

q x = {q,H 1 } D , p x = {p,H 1 } D (37) 
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and 



q x = {q, H 2 }d, Px = {p, H 2 }d, 

respectively. Here 
1 



H 1 = -((p,p) + (Aq,q)) 



and 



H 2 = ^(A 2 q,q) + ^(Ap,p) + -{Aq,q){ P ,p) -^(p,p) 2 - ^(Aq,q) 2 . 

It is easy to show that 

{A(X),A(fi)} D = {B(X),B( fJ ,)} D = 0, 

{C(A), C{ji)} D = 4(A(X) - A{p)) + 4(A^)C(X) - A(A)C(m)), 



{A(\),B(jj)} d 
{A(X),C(^)} D 
{B(X),C(h)} d 



X — fx 
2 

fJL — X 

4 

fl — X 



(B(X) - B{fx)) - 2B(X)B(fu), 

(C(X)-C(fi))+2B(X)C( l x)-2B(X), 

{A(jj.)-A(X))-4B(X)A{fi), 



(38) 
(39) 
(40) 

(41) 
(42) 

(43) 
(44) 
(45) 



where A, \x are two arbitrary parameters. 

We make use of the familiar notation [|ll]]: Li(A) = L(X) ® /, L 2 (fi) = I (8> L(fi), where 
/ is the unit 2x2 matrix; and set 



{Li(A),L(/i)}-^' mri = {L(Xy m ,m kn } D . 



(46) 



It follows from ©-© that 
Theorem 3. The Lax operator L(X) admits the following r-matrix representation 

{Li(X),L 2 (fi)} D = [r 12 (A, M ),L!(A)] - [r 21 (X, /x), L 2 ( M )] (47) 



ri2(A,/i) 
»"2i (A, At) 



At — A 
2 



At — A 



P + 2S + Q12, 
P + 25 + Q 2 i, 



(48) 
(49) 



where 



P 



( 1 \ 

10 

10 

V o o o i / 



s 





1 





1 



(50) 
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f 1 0\ / B(ji) \ / 0\ / B(u) \ 



J \ -2A(p) -B(n) J \ 1 / V ~ c (m) 

/ 5(A) \ / 1 0\ / B(X) \ ( , 
Q21= oi/u om ® n n + ^ n ® 1 n ■ ( 52 ) 



M 



2.4(A) -5(A) y V / V ~ C ( X ) / V 1 

In order to show that r%2 is an r-matrix, we must check that the bracket (24) is skew- 
symmetric and satisfies the Jacobi identity. The skew-symmetry property follows from 
the relation: 

Q 12 = -[M,L 2 (fj,)], Q 21 (X) = [M,L 1 (X)} and [[M,L 2 ],L 1 ] = [[M,L 1 ],L 2 ], 

where 

( \ 
10 

-looo • 
V o ooo/ 

Through some lengthy but straightforward calculations, we can obtain the Jacobi identity 

[Lx,[r 12 ,r 13 ] + [r 12 ,r 23 ]] + [l x , {L 2 , r 13 } d - {Z 3 , r 12 } D ] + cycle perm. = 0, (54) 

where all quantities are understood to Q <S> Q <8> Q, Q = sZ(2) (8) C(A, A -1 ); and as usual the 
superscript denotes the space which the corresponding tensor acts on nontrivially, 

Li = L <g) I <g) I, r 12 = ri2®/, etc. (55) 

An immediate consequent of (24) is that 

{L 2 (A),L 2 (^)} D = [fx 2 (X,fi),Lx (A)] - [r 21 (X^),L 2 (p)] (56) 

with 

l l 

ny(A, M ) =^^^- & (A)4-'( M )r ij (A, M )^(A)4(M). (57) 
fc=0 z=o 

Then it follows from (31) we have 

{TrL 2 (A),TrL 2 (^)} D = ^Tr {L 2 (X) , L 2 D = (58) 

which ensures the involution property of the integrals of motion obtained from expanding 
_L 2 (A) in powers of A. 

For our Lax pair L(X), we have 

N 2 N N ' , s 2 

jb=i K k=i K i=i K 1 

N ' 

Here and hereafter, . . . denotes the index of sum, which is not equal to k. It follows 

i=i 

that 

t 2 , (QkPl - Pkqi) 2 i. 1 o *r f«n\ 

J k = q k + l^ — \ x ' fe=l,2,...,JV (60) 

1 — ' Au — A] 

1=1 K 1 
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are N involution systems, which are just the conserved quantities obtained by Moser and 
Uhlenbeck j?]]. Therefore, it follows from (33) that 

Uk,Ji}D = 0, k,l = 1,2,..., N. (61) 

N 

As Jfc = (q,q) = 1, we know only N — 1 J^'s are functionally independent among 

k=l 

J\, J2, ■ ■ ■ , Jn- 
Set 

N 

3=1 

then Ft and Jj are in involution. 
Noticing 

H Q = \f x (63) 



and 



we know that both the system (|20f) and (|29|) is completely integrable. 



H t = \F 2 - l -Ff (64) 



4 Conclusion 

In this paper, we obtained the Lax representations and the r-matrix for the restricted 
KdV flows with the Neumann constraints. This approach can also be applied to other 
finite dimensional integrable systems with the Neumann constraints. Besides, for all the 
restricted KdV flows which were discussed, we can perform the separation of variables to 
integrate them following [10|; thus obtain the finite-band solution of the KdV hierarchy. 
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